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Abstract: The Gamma-Ray Large Area Space Telescope (GLAST), scheduled to be launched in Fall 
2007, is a next generation high energy gamma-ray observatory. The Large Area Telescope (LAT) instru- 
ment on-board GLAST with a wide field of view (> 2 sr), large effective area and 20 MeV to >300 GeV 
energy range, will provide excellent opportunity for future Dark Matter studies. We present an overview 
of the GLAST Dark Matter and New Physics Working Group efforts in the study of the LAT capability to 
detect a gamma-ray flux coming from WIMP pair annihilations in diffuse astrophysical sources. Particu- 
lar attention will be given to extragalactic diffuse gamma-ray radiation and line searches from annihilation 
into gamma-gamma and/or gamma-Z final states. 
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Introduction 

The nature of the Cold Dark Matter (CDM) has 
been a subject of special interest to high-energy 
physicists, astrophysicists and cosmologists for 
many years. It is probably one of the most out- 
standing open questions in present day Cosmology. 

The existence of Dark Matter is required from 
galactic scale up to cosmological scale by a wealth 
of observations and arguments such as rotation 
curve data, cluster dynamics or gravitational arcs. 
At the same time, the Cosmic Microwave Back- 
ground measurements limit the contribution from 
ordinary baryons so that non-baryonic matter has 
to make up most of the matter in the Universe. 

Virtually all potential dark matter candidates re- 
quire physics beyond the standard model of par- 
ticle physics. One of the most widely studied 
model is supersymmetric extensions of the stan- 
dard model which provides a natural candidate for 
CDM in the form of a stable Weakly Interacting 
Massive Particle (WIMP), the Neutralino. 

The mutual annihilation of these WIMPs would 
yield (among a few other indirect signatures like 
energetic neutrinos, antiprotons or positrons) many 
high energy gamma rays (>1 GeV) that can be 
well measured in the GLAST LAT. WIMP anni- 
hilation into 77 and or 7Z final states would give 



monochromatic gamma rays with an energy E~y 
m x and E 1 ~ m x (1 — m 2 z /Am x ). Since these 
gamma rays are monochromatic and have high en- 
ergy, they would constitute a spectacular signature 
of annihilating dark matter. 

The LAT Dark Matter and New Physics Working 
group has been developing approaches for the in- 
direct detection of both diffuse and point like [9] 
source observations with the GLAST-LAT detector 
[4]. We will focus here on diffuse source observa- 
tions and, as an example, particular attention will 
be given to extragalactic diffuse gamma-ray radia- 
tion [1] and line searches [6]. 



Extragalactic diffuse gamma- ray radia- 
tion 

It has been first pointed out by Bergstrom et al. 
[2] that the integrated signal of WIMP annihila- 
tions into high energy photons over all cosmolog- 
ical dark matter halos at all redshifts might con- 
tribute to the extragalactic diffuse gamma-ray radi- 
ation. 

The resulting differential energy flux per unit area- 
time-energy and solid angle on the sky can be writ- 
ten as [16]: 
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where we can distinguish contributions from parti- 
cle physics, astrophysics and cosmology. 

The particle physics contribution lies in the anni- 
hilating cross section er, the WIMP mass m x , and 
the differential gamma-ray yield per annihilation : 
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Here, dN cont / dE is the mean number of photons 
due to WIMP annihilations summed over all rele- 
vant channels and gives rise to a continuous energy 
spectrum. This term is reduced to a delta function 
in the case of monochromatic photons as expressed 
in the second term where 6 77 is the branching ratio 
into 77. 

As the annihilation rate is proportional to the dark 
matter density squared (p being the present day 
mean density), it depends strongly on how dark 
matter is distributed on small, galactic and sub- 
galactic scale. This question is still a matter of de- 
bate however N-body simulations show that large 
structures formed by the successive merging of 
small substructures, with smaller objects which are 
usually more dense [8]. The "dumpiness" of dark 
matter can then significantly boost the annihilation 
signal from cosmological WIMPs. The quantity 
A 2 describes the average of over density squared 
in haloes, as a function of redshift. 

Since we consider contributions from annihilations 
at high redshifts, the extragalactic gamma-ray sig- 
nal is strongly affected by absorption in the inter- 
galactic medium, especially at high energies. The 
term e~ T ^ z ' E °^ accounts for absorbtion of gamma- 
rays along the line of sight where the absorption 
is parameterized by the the optical depth r. We 
include the effect of absorption using parameteri- 
zations of the optical depth as function of both red- 
shift and observed energy [12]. More recent cal- 
culation of the optical depth [14] do not alter our 
results significantly. 

The last contribution comes from cosmological 
model : the Hubble parameter at the present day 
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Figure 1: 

value H , and the dimensionless quantity h(z), 
which depends on the energy content of the uni- 
verse at a given redshift. For the fractions of crit- 
ical density given by matter, vacuum energy and 
curvature, we used the results from the WMAP 
three-year data [13]. 

To compute the GLAST LAT sensitivity to Dark 
Matter contribution to the extragalactic diffuse 
gamma-ray radiation, we used the Observation 
Simulator (ObsSim) [11]. LAT photons for a 
generic model of WIMPs annihilating into bb and 
into 2 7 were generated for different WIMP masses 
ranging from 50 GeV to 250 GeV. 

Assuming that the background consists of unre- 
solved blazars [16], a % 2 analysis was performed, 
to obtain a 3 a exclusion curves in the av vs m x 
plan for one year of simulated data (fig. 1). 

Two different halo profiles have been considered 
for the normalization : the Navarro Franck White 
smooth profile [10] (NFW) and NFW profile where 
we also have included the effect of substructures. 
We assumed that substructures constitute 5% of the 
mass and have three times the concentration pa- 
rameter of the parent halo. The concentration pa- 
rameters, as a function of halo mass, is distributed 
according to [3]. 

As shown in fig. 1, depending on the Dark Matter 
halo model, the GLAST LAT would be sensitive to 
annihilating cross-sections from few 10~ 26 up to 
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One should note that this estimate neglects contri- 
butions of instrumental background, uncertainties 
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introduced by the analysis (where point-sources 
and galactic diffuse emission have to be taken into 
account) and finally that the extragalactic back- 
ground spectrum from astrophysical sources is 
very uncertain, especially at high energies. 

Line searches from WIMP annihilation 
into 77 and/or 7Z final states 

The detection of monoenergetic photons from mu- 
tual annihilation of WIMPs would be a very dis- 
tinct "smoking gun" signal of Dark Matter content 
in the Universe. These processes are however loop 
suppressed and, in popular SUSY models, these 
lines only occur at the ~ 10~ 4 to ~ 10~ 2 branch- 
ing ratio level, leading to a great challenge in de- 
tection. 

The LAT Dark Matter and New Physics Working 
group computed a preliminary estimate of the LAT 
sensitivity at 5 a above background and 95% up- 
per limits confidence level (ULCL) to these WIMP 
lines. Two cases have been investigated, the first 
being when the line energy is known (e.g. sup- 
plied by SUSY searches at the LHC). The second 
case, when energy is unknown, is less constraining 
as further statistical treatment is needed to account 
for the number of energy bins searched over the 
range of interest. 

The gamma-ray line signal depends strongly on 
Dark Matter halo and on the diffuse background 
models. The best signal-to-noise is expected [15] 
for an observation of a galactic centered broad bro- 
ken annulus (r e [25° , 35°]), which excludes the 
region within 10° from the galactic plane. 

For this study we estimated the galactic diffuse 
background in this anulus from Galprop [7] which 
is based on EGRET and other data in the EGRET 
energy range. The high energy range (up to 150 
GeV), for which the line sensitivity is computed in 
this proceeding, has been obtained from extrapola- 
tions provided by the Galprop team. 

To estimate the LAT sensitivities to known line en- 
ergy, we first considered a point source at high lat- 
itude (I = —76°, b = 26°) with a narrow Gaus- 
sian energy distribution {cf/Eq = 10~ 3 ). The In- 
strument Response Functions (IRFs) from the LAT 
Data Challenge-2 (DC-2) [5] were used, and LAT 
resolved lines were simulated with ObsSim for a 



55 day full sky scanning mode exposure (uniform 
all-sky coverage for the LAT in scanning mode). 
Energies of 25, 50, 75, 100, 125, and 150 GeV 
were considered. NB: The DC-2 version of the 
LAT IRFs cut off at 200 GeV, the currently avail- 
able LAT IRF version extends beyond 300 GeV. 

The ObsSim results for the all-sky scan were well 
fit to a double Gaussian distribution, $1 : 
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where N T = iVi + N 2 , and r = N 2 /N T . 
The all-sky diffuse background was generated for 
5 years and the flux in the galactic anulus was fit 
with an exponential background, $ 2 : 
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E/b for E e [40, 200] GeV. 



For each line, the input background was boot- 
strapped with a MC signal 1000 times and fit to 
$1 + $2 and $ 2 , to compute the 5a signal sen- 
sitivity for 5 years of observations. This series of 
1000 bootstraps was re-run varying only the aver- 
age number of MC signal counts (constant across 
the series) until < Ax 2 >= x! 1+ <e> 2 - X% 2 ~ 25 
(corresponding to a 5 a signal over background). 
The average number of signal counts needed at 
each energy was then converted to a sensitivity us- 
ing average exposures integrated over the annulus 
(fig 2). The LAT detection sensitivity for a line 
of unknown energy (E € [40, 150] GeV) for a 5 
a above background signal corresponds to a con- 
fidence level of 10~ 7 . To calculate the number of 
counts needed at each energy in this case we used 
the probability of no detection in a single bin, q: 

q = (l- p)V#ofEnergyBins 

= -±= / e-* a /2 dx, 

where P = 10~ 7 is the probability of detecting a 
signal greater than #o/cr, in one or more bins. A 
bin width of <Jf w hm/E = 8% was used based on 
the FWHM energy resolution of the IRFs. The er- 
ror in the signal was estimated as a S i gna i = V2B 
where B is the number of background counts in a 
bin centered on Eq with width (Jf w hm- The #o/ct 
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Figure 2: 

signal error was converted to flux for calculation of 
the unknown line energy (fig. 2). 

The 5 year 95% CLUL sensitivity for known line 
energies was obtained similarly to the 5a case by 
bootstrapping the ObsSim diffuse background and 
fitting $i + $2 with no MC generated signal. Af- 
ter the calculation of N T ± <jnt, the 95% CLUL 
was obtained from 1.64 <jnt and sensitivities cal- 
culated as in the 5er case with P=.05 (fig 3). The 
case of unknown line energy (E 6 [40, 150] GeV) 
95% CLUL was calculated as in the 5a case with 
P = .05 ( fig 3). 

Conclusions 

We presented the GLAST Dark Matter and New 
Physics Working Group efforts on indirect detec- 
tion of Dark Matter through diffuse source obser- 
vations. As an example, preliminary results on 
extragalactic diffuse gamma-ray radiation and line 
searches illustrated that the GLAST-LAT detector 
will provide excellent high energy gamma-ray ob- 
servations for future Dark Matter searches. 
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